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Abstract 



We study fermionic excitations in a cold ultrarelativistic plasma. We con- 
struct explicitly the quantum states associated with the two branches which 
develop in the spectrum as the chemical potential is raised. The collective 
nature of the long wavelength excitations is cleary exhibited. 
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1 Introduction 



The spectrum of fermion excitations in a hot relativistic plasma composed of elec- 
trons, positrons and photons (or quarks, antiquarks and gluons) presents an inter- 
esting structure [1-8]. In particular, for small fermion masses m and momenta p, 
i.e. m,p ^ gT where g is the coupling constant and T the temperature, it is split. 
One of the two branches becomes the ordinary fermionic branch at large momentum 
or large mass. The other has features of a collective excitation, and for this rea- 
son has been sometimes dubbed the "plasmino" B, by analogy with the plasmon 
oscillation. However, one should keep in mind that, when p < gT, both branches 
have a collective character and, as we shall see, the whole structure of the spectrum 
at small momentum reflects a collective behavior of the system. A recent calcu- 
lation indicates, at least within the one-loop approximation, that this structure in 
the spectrum does not occur as an abrupt transition, but develops gradually as the 
temperature is raised ||. 

Similar effects are expected to occur at zero temperature, but large chemical 
potential //, i.e. when g\i ^ m,p@. This is the situation we explore in this paper. 
There are many physics motivations for studying the spectrum of fermions in very 
dense matter, ranging from the physics of dense stars, to that of heavy ion colli- 
sions. These have stimulated recent worksJTDJ. There are also current investigations 



on more formal questions related to QED with a chemical potent ialJTT]]. However, 
the main purpose of the present paper is simply to get a better physical under- 
standing of the fermion spectrum, especially at low momentum where collective 
phenomena develop. We do this, in particular, by providing an explicit calculation 
of the quantum states corresponding to the fermionic excitations, which is possible 
at zero temperature. 

Some time ago, Weldon presented a physical interpretation for the splitting of 
the spectrum at high temperature, based on an analysis of the quantum numbers of 
the modes[]l2|], and an analogy with BCS theory of superconductivity ||] . However, 
that did not totally clarify the simple physics underlying the emergence of collective 
fermionic excitations in the high temperature plasma. In fact the analogy with BCS 
superconductivity is misleading. The formation of Cooper pairs involve particles 
localized in momentum space at the vicinity of the Fermi surface, and it produces 
an instability of the Fermi sea leading to a reorganization of the system as a whole. 
Here, as we shall see, the particles in the vicinity of the Fermi surface are only 
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weakly perturbed, and the Fermi sea is stable. 



We begin in the next section by summarizing some properties of the fermion 
self-energy at finite temperature and chemical potential, to leading order in the 
coupling strength g. The spectrum of fermion excitations is presented, and some 
known results are recalled. We show in particular that, in the ultrarelativistic limit, 
most properties of the spectrum depend on [i and T only through the combination 
M 2 oc <? 2 (/z 2 + 7r 2 T 2 ). In section 3, we discuss some general properties of the one- 
loop fermion self-energy at zero temperature and finite chemical potential. Analytic 
results are obtained in two simple cases, m = and p = 0. We discuss the approach 
to the ultrarelativistic limit and the emergence of a particle-antiparticle symmetry 
at large p. We show that, as it is the case at finite temperature, the structure in 
the spectrum of quasiparticles develops gradually as the chemical potential is raised, 
due to the coupling of the quasiparticles to low lying states involving hard electrons 
and photons. We argue that the imaginary part of the self energy may not be 
interpreted as a damping of the quasiparticles. The main results of this paper are 
contained in section 4 where we present a detailed calculation of the quantum states 
corresponding to the two branches in the excitation spectrum at small momentum. 
The collective nature of the two branches is clearly exhibited, and various features 
of the spectrum are analyzed. Our conclusions are summarized in the last section 
which contains also a discussion of various aspects of the present work, such as 
its relevance to the case of finite temperature, and also to QCD, the relation with 



kinetic theory developed in [13|], and finally the relation with a mode analogous to 
the plasmino, the "plasmaron" found in some calculations of the non relativistic 
electron gasO, ITH [16 . 



2 Fermionic modes in ultrarelativistic plasmas 

In the presence of a chemical potential associated with the conserved charge Q = 
J ipjoi/jcfix, the Dirac hamiltonian en-p + rrryo becomes a -p + 77170 — p ( ex = 707, 70 
and 7 being the usual Dirac matrices). Thus, the net effect of the chemical potential 
is to shift all the single particle energies by — p. Correspondingly, the inverse of the 
free Dirac propagator may be written 

Go X (p,uj) = -(^ + /i)7o + 7-p + m 

= (e P - - ^hoAp - (jEp + p + w)7oA p (1) 
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where 



A p = 2^ (ep±( "' P + m7o)) (2) 

are the projectors on positive and negative energy solutions respectively ((A 1 * 1 ) 2 = 
A ± , A + + A - = 1), e p = y/p 2 + rri 2 , and uj is a complex variable whose real part 
measures a single particle energy with respect to the chemical potential. In the 
following, we shall in fact absorb \i into a>, i.e. we substitute uj to uj + \i in Eq.fll]) so 
that Re uj will measure the energy with respect to the zero of energy in the absence 
of chemical potential. If needed, the Feynman propagator is obtained as usual by 
letting uj approach the real axis from above when Re uj > //, and from below when 
Reu </i§|. 

For simplicity, we consider in the main text the interaction of a Dirac fermion 
if) and a massless scalar field (ft with a Yukawa coupling g<jnpip. However, most of 
the discussion extends to QED (see Appendix B) and for this reason we shall often 
refer to the fermion as an "electron" and to the scalar particle as a "photon". In 
presence of the interaction, the inverse propagator is G^ 1 = Gq 1 + E, and, owing to 
rotational invariance, the mass operator E may be written in the form 

E(p, uj) = a(p, uj)<y° + b(p, u)*y ■ p + c(p, uj) (3) 

where a, b and c are three independent functions of the length p of the momentum, 
and the energy uj. 

To order g 2 the mass operator at finite temperature and chemical potential is 
the sum of two terms E(p, uj) = Ei(p, uj) + E (p,c<j), where E (p,co | ) is the vacuum 
contribution, that is the mass operator for T = \i = 0, and E^p, uj) is the correction 
due to finite /i, finite T, effects. The latter is given by 



where 



o / d 3 k 1 / A+ n k - / p+k A _ /p+k + n k 
(2n) 3 2uj k \ p+k cj-e p+k -^ p+k cj + e p+k - u k 

+ A p+k — - — - — + A P +k— — no, (4) 

uj - e p+k + uj k uj + e p+k + u k .' 



e ^(6 P ±M) + i nfc e ^ - 1 ^ 
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are the occupation numbers respectively for the electrons (fp), the positrons (fp), 
and the photons (n&). The photon energy is u>k = k, except in the calculation of 
the vacuum contribution (see below and Appendix A) where we set uj k = Vk 2 + A 2 , 
A playing the role of an infrared cut-off. Note that in writing Eq.(f|), we ignore 
contributions involving possibly non-vanishing (in fact infrared divergent) expecta- 
tion values of iptp (or in the case of QED). However, such contributions to S 
are independent of energy and momentum and produce simply a constant shift in 
the single particle energies. Therefore, they can be absorbed in a redefinition of the 
chemical potential. Now, since this plays no role in the present discussion, we simply 
ignore these constant contributions. (In the usual treatment of the non-relativistic 
electron gas, one assumes that such contributions are cancelled by those arising from 
a positive background of positive charges. In the case of QCD, such contributions 
vanish because of color neutrality). Thus, in zeroth order, the chemical potential 

will be taken equal to the Fermi energy, i.e. /x = \jm 2 + hp, where k F denotes the 
Fermi momentum, related as usual to the density of the conserved charge. 

The vacuum (T = /x = 0) contribution to the mass operator has a structure 
similar to Eq.(f§), 

v ( \ 2 f d3k 1 f A + 1 a- 1 \ 

S ° (P)U;) " 9 J (^F^i Ap+k c - e p+k - u k + + e p+k + u k P°- 

(6) 

The integrations over momentum are diverging and require renormalization (see 
Appendix A). However, as we shall see shortly, in the ultrarelativistic limit where 
m ~ g/j,,gT, and for soft modes with u>,p ~ gfi,gT, the contribution of S is of 
higher order in g than the dominant terms of S 1; and it will not be discussed further 
in this section. 

In the limit where m, u, p /x, T, the integrand in Eq. (f|) contributes mostly 
in the region k ~ /x or T where we can replace e p+ k by k +pcos8, with 9 the angle 
between p and k. Furthermore, the first and the last term in Eq. (f|), for which the 
energy denominators are of order k ~ /x, T, are small compared to the second and 
third terms where the energy denominators are approximately uj ± p cos 9 T, /x. 
The same reasoning shows that the vacuum contribution, which has the same denom- 
inators as the first and last terms of Si, can be discarded in a first approximation. 
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Keeping only the dominant terms, one obtains then 



E(p,w) 



d 3 k 1 



A 



P+k" 



f + 



UJ k 



16tt 



(2ixf2uj k l uJ + e p+k 
^ MA; + f+ + 2n k ) j ^ 



+ A P+k 

g?cos6* 



/p+k + 



e p+k + . 



UJ 



pcos9 



1 



1 ■ P 

P 



cos 9 



(7) 



The integrals over k and over 9 are decoupled and can be done exactly. Note that the 
fc-integration extends to the region of soft momenta for which the approximation 
done is non longer valid. Therefore, the expression (|?p should be used only in 
situations where the weight of such soft momenta is negligible. A straightforward 
calculation gives the following results for the functions a, b and c defined in Eq.(RI): 



with uj real and 



Re a(p, uj) 


M\ 
= In 

2p 


uj + p 
uj — p 


Re 6(p, uj) 


p V 


UJ 

~ 2p 


Re c(p, a;) 


= 




M 


= A^ 2 + 7r2T2 . 



In 



UJ + p 



uj — p 



(9) 



In the ultrarelativistic limit, the self-energy E depends on the chemical potential 
and the temperature only through this parameter M. 

The functions a and h have imaginary parts when \uj\ < p, given by 



Im a(p, uj) 
Im b(p, uj) 



— 7T- 



~2p~ 

M 2 uj 



71- 



2p 3 



(10) 



if — p < uj < p and otherwise. The states to which the soft fermion couples, 
and which are responsible for this imaginary part, may be associated to virtual 
transitions where hard particles scatter on the soft one, with little deflection, and 
a possible change in their quantum numbers. For example, a hard electron can 
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annihilate on a soft positron and "turn into" a hard photon. The kinematics of 
such processes is somewhat analogous to that involved in the Bremsstrahlung of 
soft photons by fast electrons. Here, the soft particles are electrons or positrons 
and the hard particles are electrons, positrons or photons. These virtual transitions 
have uj ~ pcosO and their number is constant, independent of p; thus their density 
increases asl/p when p decreases. In fact, this density is measured by the imaginary 
part of a oc M 2 /2p. Essentially all particles can contribute to these transitions, as 
reflected in the momentum integral in Eq. (|7|). Note that when p — > 0, b — > 0, 
and a(p,u) — > M 2 /uj, so that, in this limit, Ima(w) oc 5(uj). The existence of an 
imaginary part in the region \uj\ < p is reminiscent of the familiar Landau damping 
of electromagnetic waves in ordinary plasmas. Here, however, the imaginary part is 
of the same order of magnitude as the real part, i.e. O(M), and excitations with 
such energies (\uj\ < p) would be quickly damped. 

Well defined quasiparticles exist for \u\ > p. Their energies are the poles of the 
propagator, and are therefore given by the solutions of the equation det(G _1 ) = 0, 
where 

G-\p, u) = (a - u) 70 + (1 + b) j ■ p + (c + m), (11) 

or 

uj - a = e + bf p 2 + (c + m) 2 , (12) 

with e = ±1. Using the explicit form of the functions a, b and c just derived, one 
can write Eq.(^) in the form 




with 

x + 1 
h(x) = 1 — In 

The solution of Eq. ([13]) gives Mjp as a function of uj/p, from which one obtains 
the dispersion relation uj(p) in parametric form for a given M. Note that Eq.(|l3D 
is even in uj, and its solutions come in pairs with opposite signs. Thus, in the 
ultrarelativistic limit, there is an apparent particle-antiparticle symmetry for the 
soft quasiparticles excitation energies. This is at first sight surprising since such a 



x+ 1 



X 



(14) 
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symmetry is explicitly broken by the presence of the chemical potential. We discuss 
this point further in the next section. We note here that if u is a solution of Eq. (|13| ) 
with e = 1, then — u is a solution with e = —1. Furthermore, for a given e, there are 
two solutions of opposite signs. We call uj + and u;_ the two solutions corresponding 
to e = 1, with uo + > 0, lo_ < 0. We call "normal" the branches eu + and "abnormal" 
the branches £a>_. Thus, for u > 0, the branch u + is the normal one, the branch 
— co>_ the abnormal one. This seemingly complicated notation will be fully justified 
in section 4 after we have clarified the origin of the observed pattern of states. As 
we shall see, e can be identified with the eigenvalue of 70 in the case where p = 0, 
and with that of a ■ p when m = (p = p/p). 

The spectrum is displayed in Figs.l. If m is small enough compared to M, the 
group velocity on the abnormal positive branch, — du^/dp, is negative for small p, 
and — uj- (p) presents a minimum for a finite value of the momentum. The condition 



under which this phenomenon occurs, m/M < yl4/(3\/5) — 2 w 0.295, can be 
derived easily by expanding Eq . (|i~3"D in powers ofp/u. For m = 0, we have 



u±(p) = ±M 



l±i^ + ... 
3M 



(15) 



Thus, in particular, the slopes of the two branches uj± are the same at small p, 
duj±/dp = 1/3 when m — 0. As M becomes smaller than m, the normal branch 
eu + (p) goes over to the normal fermion dispersion relation Ey/p' 2 + m 2 , while the 
abnormal branch approaches — ep. The asymptotic behavior forp ^> M and 

m = 0, is 

M 2 

w+(p) ~ p+ 

P 

uj-(p) ~ -p [l + exp(-2p 2 / M2 )l (16) 

The residues at the quasiparticle poles can also be calculated. They have simple 
expressions in the two limiting cases where p = or m = 0. Then, the projectors 
(0) reduce to (1 ± 7o)/2 when p = 0, and to (1 ± a. • p)/2 when m = 0, and we can 
write , in the appropriate subspace, 

G(p, ^70 = ^— — ss — — — (17) 

{ee p — uj) + l £ {uj) euj± — uj 

where e denotes the eigenvalue, ±1, of 70 in the case p = 0, or that of ct ■ p when 
m = 0, S e = a + ec when p = 0, and S e = a + ebp when m = 0. The residues are 
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given by z± = 1 — dY* £ / du\ UJ± and they depend only on the character of the mode, 
normal or abnormal, irrespective of the sign of u, i.e. of e. When m = 0, 



When p = 0, 



2± = T7W (19) 



The variations with p of the residues are displayed in Figs.l. When m 7^ 0, the 
residue of the abnormal branch is always smaller than that of the normal branch. It 
becomes negligible when either m > M or p > M. For p = we have z + + z_ = 1. 
For p^O, z + + z_ < 1; this is because the spectral weight has a smooth component 
between uj = —p and u = p. We will comment on this in section 4. 



3 One loop results at T = 



We now study the full one loop correction to the fermion propagator at T = and 
p 7^ 0. In this case, the mass operator has a simple analytic expression when either 
p = or m = 0. The corresponding formulae are collected in Appendix B. Some 
of the discussion in this section parallels that in Ref.||, where a similar study was 
presented at finite temperature and zero chemical potential. It is instructive to see 
how the structure in the spectrum develops as the chemical potential is raised and 
the ultrarelativistic limit is approached. However, one should keep in mind that the 
full one-loop calculation is not entirely consistent, as revealed for example by the 
order of magnitude of some contributions to the imaginary part. We shall discuss the 
physical meaning of these imaginary parts and see why it is not always meaningful 
to associate those in the vicinity of the quasiparticle poles with the damping of 
quasiparticles. 



3.1 The mass operator at T = 

At zero temperature, there are no bosons, i.e. = 0, nor antifermions, i.e. = 0, 
and all single particle states with energies less than p, are occupied, i.e. f~ = 
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6(fi — e p ). The mass operator takes the form 

E{p,u)=g I 77^ — 1 A p+k 



(2vr) 3 2^ fc l p+k cu-e p+k -cj fc 



+ A P +k +A p+k 7 1 7 k, 

uj - e p+k + u k uj + e p+k + u k J 

(20) 

where we have kept together the vacuum contribution and the finite chemical po- 
tential correction, without writing explicitly the renormalization counterterms (see 
Appendix A for details). Note that the imaginary part, which we discuss below, is 
finite and not affected by the renormalization. In order to facilitate the foregoing 
discussion, we label by (i), (ii) and (iii) the three contributions to E(p, uj) in Eq. (p0|). 
These are illustrated by the diagrams of Figure 2. In these diagrams, the labels on 
the lines correspond to the momentum and the energy of the single particle state to 
which the line refers. A fermionic line going upwards represents an electron above 
the Fermi surface. A fermionic line going downwards represents a hole in the Fermi 
sea if the energy is positive (case (ii)), and a positron (a hole in the Dirac sea) if the 
energy is negative (case (iii)). Thus, reading the diagram (ii) in Fig. 2 from bottom 
to top, one finds a hole excitation with energy — uj and momentum — p which makes 
a virtual transition to a state containing a photon with momentum k and energy k, 
and a hole in the single particle state with momentum p + k and energy e p+k . 

It is perhaps useful at this point to recall that the poles of the Green's function 
occur for values of uj equal to ±(E^ ±1 — Eq), where E^ 1 is the energy of a 
state with N ± 1 particles, and total momentum p, while Eq is the energy of the 
ground state with N particles. (Note that in this discussion we count particles with 
respect to the fully occupied Dirac sea, and similarly for the energies.) Perturbative 
calculations are meaningful as long as Eq — Eq' 1 < fi < Eq +1 — Eq fl8fl . Values of 
L) such that Re to > fi correspond to possible energies of states with + 1 particles, 
while values of uj such that Re to < fi correspond to the negative of energies of states 
with N — 1 particles. For example, the state obtained by removing a particle of 
momentum p in the unperturbed Fermi sea has a momentum — p, and an energy 
— uj = E^ 1 — Eq = — e p , with uj < ft. The same reasoning applies to positrons 
considered as holes in the Dirac sea. Removing an electron in the Dirac sea in a 
single particle state of momentum p and energy — e p , produces a state with — 1 
particles, a momentum — p, and an energy — (— e p ) = e p . 
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The mass operator acquires imaginary parts in several energy domains where the 
energy denominators in Eq. (|20| ) vanish for uj real. To discuss these, we note the 
following inequalities 

e P +k + k > e p , < e p+k - k < e p . (21) 

It is then easily seen that, quite generally, the imaginary part of (i) vanishes when 
uj < /i. This follows from the fact that the occupation factor in the numerator 
forces e p+k > /i, so that the denominator uj — e p+k — k cannot vanish if uj < [A. 
Similarly, the occupation factor / p+k in the numerator of (ii) forces e p+k < p so 
that the imaginary part of (ii) vanishes when uj > fi. Finally, using the first of the 
inequalities (|2T| ) one gets that the imaginary part of (iii) vanishes when uj > —m. 

Now, the precise domains in which the imaginary part does not vanish depend 
on p. These are illustrated in Figure 3. The last term of Eq. fl20f) , (iii), contributes 
an imaginary part when uj < — e p , as easily deduced by using the first of inequalities 
(pTT) ■ This corresponds physically to the process by which a positron (off its energy 
shell) with energy \u\ and momentum p decays into a positron of momentum p + k 
and a photon of momentum — k, on their respective mass shells; this process may 
occur in the vacuum, for example under the action of a weak perturbation which 
displaces the positron slightly off its mass shell, and it is not affected by the Fermi 
sea of electrons. We shall refer to the continuum of states into which a positron can 
decay as the positron-photon continuum. It corresponds to region (III) in Fig. 3. 

The first term in Eq. (|20D , labelled (i), has, for p > kp, an imaginary part when 
uj > e p . This imaginary part lies in the vicinity of the electron mass shell, and reflects 
the possible decay of an electron above the Fermi surface into another electron with 
less energy and a photon. As in the positron case just discussed, this process may 
also occur in the vacuum; however, here, the presence of the Fermi sea eliminates 
some of the states available for the decay in free space, because of the Pauli principle 
reflected by the occupation factor 1 — / p+k in Eq. (^Cj). In the case p < k F , this 
factor actually reduces the energy domain allowed by the kinematics, i.e. uj > e p , 
toco> > fi + kp — p, and indeed leads to a cancellation of the imaginary part in the 
region e p < uj < fi + kp — p. Then, for p < kp, the imaginary part of (i) is non zero 
in a domain (u > fi + kp — p) which can be far from e p and it is not immediately 
related to a possible decay of the single particle excitation. The continuum of states 
discussed in this paragraph will be referred to as the electron-photon continuum. It 
corresponds to region (I) in Fig.3. 
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The existence of states of the N + l particle system with momentum less than the 
Fermi momentum kp reflects modifications of the Fermi sea due to the interactions. 
Because of these, the single particle states are neither fully occupied below the 
Fermi surface, nor completely empty above it. These modifications of the Fermi 
sea are implicitly taken into account, to order g 2 , in the expression ( ^0[) of the mass 
operator. The interplay between the single particle properties and these ground state 
corrections is perhaps best understood by considering the energy shift of a single 
particle state obtained in second order perturbation theory. This is obtained by 
solving Dyson's equation for the pole of the propagator, setting u = e p in the mass 
operator E(p, uj). The first term in S(p, e p ), i.e. (i), gives an attractive contribution, 
while the second term, i.e. (ii), gives a repulsive contribution. This is easily seen 
by using inequalities (0): the denominator of (i) is always negative, while that of 

(ii) is always positive. To understand this better, consider a particle state above the 
Fermi surface, i.e. with p > kp, and recall that the poles of the Green's function 
correspond, for ui > /i, to the energy difference E^ +1 — Eq, where E^ +1 is the energy 
of the system with N + 1 particles and momentum p. In carrying out the second 
order calculation of this energy difference, one is lead to write the total correction 
as Ae p + AEq +1 — AEq . The first term, Ae p , is the correction obtained from (i); 
it is negative and can indeed be interpreted as the second order correction to the 
energy of the single particle state p, due to its coupling to states with one electron 
(above the Fermi surface) and one photon. On the other hand, the quantity AEq +1 
represents a correction to the the energy of the Fermi sea, taking into account the 
presence of the extra particle. The full second order correction is given by AEq , 
and is negative. Because of the Pauli principle, the presence of the extra particle 
suppresses some of the intermediate states (those having a particle with momentum 
p above the Fermi surface) involved in the full second order correction to E^; this 
is illustrated in Fig. 4. As a result, AEq +1 > AEq , and the contribution of (ii), 
which can be identified to AE$ +1 — AEq , is positive. A similar interpretation may 
be given for states below the Fermi surface, exchanging the role of contributions (i) 
and (ii). 

We come now to the second term in Eq. fl2"0|), labelled (ii). In contrast to (i) and 

(iii) where vacuum contributions explicitly enter, this term has no counterpart in the 
vacuum and is entirely due to the Fermi sea of electrons. It contains the dominant 
contribution in the ultrarelativistic limit. It has for p < kp an imaginary part in the 
region fi — kp—p<u)<e p . This imaginary part measures the density of states with 
momentum p containing one hole and one photon with respect to the unperturbed 
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Fermi sea. Such states may couple with a hole at the bottom of the Fermi sea, or 
a hole at the top of the Dirac sea, i.e. with a positron. A pictorial representation 
of all these states, which will play a central role in our discussion (in particular in 
section 4), is given in Fig 5. One may view the corresponding virtual transitions as 
processes by which a hole in the Fermi sea, or a positron, gets "filled up" by electrons 
falling into it from above in the Fermi sea, emitting at the same time a photon. If 
p > k F , (ii) has an imaginary part in the domain jjl — p — k F < uj < fi — p + k F . The 
existence of states of the N — 1 particle system having a momentum greater than 
the Fermi momentum reflects, similarly to the case discussed above, a modification 
of the Fermi sea due to interactions. The continuum of states corresponding to the 
virtual transitions described in this paragraph will be refered to as the hole-photon 
continuum. This correspond to region (II) of Fig. 3. 

In the ultrarelativistic limit, the hole-photon continuum is entirely contained 
in the region \u\ > p (see section 2, the discussion after Eq.flTDp. However, the 
above discussion shows that the imaginary part does not strictly vanish outside this 
interval. The additional contribution, in the interval p < uj < e p , can be calculated 
assuming \i ^> uj,p,m and uj — p g 2 ji. The result is, 

q 2 uj \m 4 — (uj 2 — p 2 ) 2 ] 
Im a = -7T- 1 ' 



Im b = — 7r 



327T 2 (uj 2 -p 2 ) 2 
g 2 m 4 — (uj 2 — p 2 ) 2 

32^2 (^2 _ p 2)2 



Imc = -*4- 2 2m[m2 ;^- p2)] (22) 
Slix 1 uj 1 — p z 

The same expressions hold also, up to a global change of sign of Imfr, in the region 
uj < — e p correponding to the positron-photon continuum. When p, m and uj are of 
order g/i, Im E(p, uj), as given by Eq. (|22"D , is of order g 3 fi, i.e. smaller by two powers 
of g than when \uj\ < p. This comes from the fact that almost all the states in the 
Fermi sea contribute to the hole-photon continuum, while only a restricted phase 
space is involved in the calculation leading to Eqs.fl2"2|). Indeed the phase space 
vanishes as uj approaches e p , as may be seen on these expressions (|2"2"D. Note that 
Eqs.(^) do not involve the chemical potential. This is to be contrasted to the case 
T ^ 0, where the phase space depends explicitly on the temperature through the 
statistical factors. Thus at T ^ the imaginary part of the one-loop mass operator 
is of order g 2 T (instead of g 3 /i here), because the number of soft photons of energy 
~ gT is of order l/g. 
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We have seen in section 2 (see also Fig.l) that, in the ultrarelativistic limit, there 
are four quasiparticles of a given momentum p <kp, whose energies are the solutions 
of Eq. (|T2"D . Some of these quasiparticles appear in regions where ImE is small but 
non vanishing. One may then be tempted to associate these small imaginary parts to 
a weak damping of the quasiparticles. Because of them, the quasiparticle energies 
are indeed shifted to complex values, to — uj± — i r y±, where 7± is always positive 
(7± = -z±ImS(u±)). (It is interesting in this respect to notice that without the 
vacuum contribution, there would be some imaginary part for e p < uo < fi coming 
from the term labelled (i), with the wrong sign.) The normal quasiparticle branch 
with e = 1 always satisfies uj + > e p and thus has no imaginary part (the branch 
uj + is above the hole-photon continuum). On the other hand, the normal branch 
— u + < — e p is always damped (it is always in the positron-photon continuum). 
As to the abnormal branches, there are two possibilities. Either — u;_ > e p (case 
M > m in Fig.l); then the branch — u_ is above the hole-photon continuum, 
and the branch uo_ < is inside the positron-photon continuum. Or — u;_ < e p 
(case M ~ m in Fig.l), and then it is the mode — uj- > which is damped, 
being found inside the hole-photon continuum. Now one may question whether 
the fact that a branch of the spectrum lies inside a continuum of states implies 
necessarily a damping of the corresponding single particle excitation. When M is 
sufficiently large compared to m, possible damping would occur because of a mixing 
of the single particle excitation with the positron-photon continuum . However, a 
consistent calculation of the damping should take into account the fact that the 
edge of this continuum, taken here to be the unperturbed positron energy e p , is 
strongly modified by the interactions, which affect actually both the propagation 
of the fermions and that of the photons. Taking as a minimal approximation the 
modification of the fermion energy in the calculation of the damping, one would 
obtain that the damping is forbiden by the kinematics (the branch — u + becomes in 
this approximation the edge of the continuum). The difficulty encountered here is 
identical to that noted in a related context (the damping of gluonic excitations in 
a quark-gluon plasma) by several authors |fTP|, |2H| , |2"T| |. Technically, it finds its origin 
in the fact that the imaginary parts, as obtained from the one loop approximation, 
being of order g 3 /i, are inconsistent since higher loop diagrams can contribute to the 



same order. Consistent calculations require resummations as explained in Refs. |22 

In order to get further insight into the properties of the quasiparticles, and the 
characteristic features of the ultrarelativistic limit, we now briefly consider succes- 
sively the particular cases where p = and m = 0. 
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3.2 Massive fermion at rest 



The equation which determines the poles of the propagator (in the region where the 
imaginary part is small), namely Eq . (|T2|) , becomes, when p = 



uj — em = Re(di + eci), (23) 

where we have ignored the vacuum contribution, and e denotes the eigenvalue of 
7°. Analytical expressions for Reai and Reci, as well as for the corresponding 
imaginary parts can be found in Appendix B. The quantity Re(ai + eci) is plotted 
in Figs. 6, as a function ofu/m for various values of p/m. In general, when /i ^ 0, 
there is no symmetry under the transformation u — ► —uj, e — > — e. However, this 
symmetry emerges at large p, as is already apparent in Figs. 6 (the function Reai 
becomes odd in u and Reci becomes negligible at large p). A graphical solution 
of Eq. ( P3|) is obtained by taking the intersection of the curve Re(ai + ec\) with 
the straight line u — em. One sees that the normal unperturbed pole at uj = em is 
shifted to the right (e = 1) or to the left (e = —1). Furthermore, when e = 1, two 
additional poles appear on both sides of the logarithmic divergence at u — p — kp, 
even if p is arbitrarily small. When e — —1, a similar situation occurs with two new 
poles slightly above to = p — kp, but only if p is large enough. The physical relevance 
of these extra solutions may be asserted by considering the imaginary part of a + ec 
(see Eq. ( [B.3| )). One sees that the imaginary part remains large in the vicinity of 
the pole closest and above p — kp. This pole correspond to an excitation which is 
killed by Landau damping. The other two poles correspond to physical excitations 
with negligible widths. As we have seen earlier, it is consistent to ignore the small 
imaginary parts in the ultrarelativistic limit. A consistent calculation of the damping 
rates of the quasiparticles requires going beyond the one loop apparoximation. 

The structure of the self-energy can be understood from the dispersion relation 

ReS w = P / ^ 24 

v ' J 2tt u'-u K J 

where the symbol P in front of the integral indicates that the principal value is to 
be taken. The region contributing to the dominant imaginary part is the small u 
region where the imaginary part goes as uj~ 3 . This behavior reflects the increase of 
the phase space for processes in which a hole at the bottom of the Fermi sea couples 
with holes near the Fermi surface and photons with momentum ~ kp. When /i 
grows, the lower bound of the allowed phase space p, — kp ~ m 2 /(2p) — > 0, and 
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the number of such transitions becomes increasingly important. Thus, at large /x, 
the imaginary part is approximately given by ImS ~ — (g 2 m 4 /327rc<j 3 ). Using the 
dispersion relation (JM|), and keeping only the dominant contribution coming, when 
uj ~ gfi, from the lower bound of the integral ~ m 2 / (2/i), one gets Re £ ~ M 2 /uj = 
Rea(0,c<j) (M = gfi/4ii), in agreement with Eq.(H). Thus the leading behavior of 
Re £ is determined by the hole-photon continuum, to which all electrons in the 
Fermi sea can contribute. This explains the emergence of the symmetry uj — > —uj in 
the ultrarelativistic limit. Most of the particles involved in the virtual transitions 
making the hole-photon continuum are hard particles. In pictorial terms, an electron 
with momentum p ~ fi ^> m falling into a hole with momentum p ~ gfi, does not 
"see" the mass gap, i.e. it does not feel the difference between a hole at the bottom 
of the Fermi sea and a hole at the top of the Dirac sea. 

In the high density limit Eq. ( p.2|) has only two solutions for either value of e 
(the third solution, slightly above fi — k F , disappears in this limit), that is to say 
four solutions. The positive solutions are 

m / „ f m\ 2 
±uj ± = ±- + ^M 2 +(-) . (25) 

Eq. ( |B.4|) shows that the imaginary part vanishes for both positive solutions if M > 
m\/2, and only for the solution e — 1 if M < mi/2. As to the negative solutions, 
— uj has no imaginary part when u has one and vice-versa. As already mentioned, 
this imaginary part is small and can be ignored. 

3.3 Massless fermion with a finite momentum 

When m = 0, Eq.(|T2"D becomes (as in the previous case we ignore the imaginary 
parts and the vacuum contribution) 

uj — ep = Re (ai + ship) (26) 

and e denotes the eigenvalue of a ■ p (p = p/p)- The exact expressions for the real 
and the imaginary parts of a\ and b\p are given in appendix B, and Re(ai + ebip) is 
displayed in Figs. 7 as a function of uj. One sees that the imaginary part of a + ebp 
is continuous everywhere except at uj = —ep, as can be checked in Eq. ( B.7|) . This 
discontinuity at the edge of the hole-photon continuum is an effect of the zero 
mass. The imaginary part is significant only if \u\ < p. The number of processes 
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contributing to it is an increasing function of \i. We also note that, as in the previous 
section, the real part of a\ + eb\p is, in general, not odd under the changes of sign 
of both e and uo. However, it becomes so in the high /i limit, for reasons which have 
already been explained. Eq. ( p6|) can be solved graphically in Fig. 7, exactly in the 
same way as Eq. (|3|) in Fig. 6. One sees that for a small value of g, Eq. (p6|) always 
has two solutions in addition to the solution close to cu — ep. These new poles are 
on both sides of u = —ep, and the pole which is inside the hole-photon continuum 
is Landau damped if /i ^> p. Thus we have finally one new physical pole for either 
value of e when m = 0, as we had for p = 0. 

Finally, let us mention that the states in the vicinity of the Fermi surface are not 
much affected by the interactions in this order. This is illustrated in Figs. 8. The 
scale on these figures shows indeed that £ is a small quantity when p ~ hp (it is 
completely perturbative, i.e. of order g 2 ). The lines w^fp are almost vertical lines in 
these plots, indicating that single particle states are shifted from their unperturbed 
positions by negligible amounts. As for the imaginary part, it vanishes for u = [i as 
is evident on Eq. (|B.7D . 

We end this section by summarizing the physical picture which emerges from this 
analysis. As we have seen, most of the interesting physics is due to the hole-photon 
continuum. When fi is large, this continuum involves virtual transitions in which 
hard electrons "fall" into holes at the bottom of the Fermi sea, or at the top of the 
Dirac sea (soft positrons). The number of such transitions grows as /i 2 , reflecting 
the fact that most particles of the Fermi sea participate and suggesting the collective 
character of the quasiparticles. We shall now complete this picture by presenting an 
explicit calculation of the quantum states corresponding to the quasiparticles. 

4 Explicit calculation of the states 

In this section we present an explicit construction of the quantum states associated 
with each branch of the excitation spectrum. We shall work mainly in the ultrarel- 
ativistic limit. As we shall see, many of the features which were left unexplained in 
the previous analysis will find here a natural interpretation. 

We decompose the fermion and scalar fields into plane waves normalized in a 
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volume f2 



= ^E^ p - x ^A,p&A, P (27) 

«*> = ^^^(e-^ + e— 4) (28) 

In the decomposition of the fermion field if}{x.), we do not separate at this stage 
positive and negative energy solutions of the Dirac equation, and consider a positron 
as a hole in the sea of negative energy electrons. In Eq.(^7]), we therefore sum over 
positive and negative energy solutions and the index A takes four values, two for the 
spin and two for the sign of energy. This convention allows us to treat the Dirac and 
Fermi seas on the same footing: all single particle states with energy less than /z are 
occupied in the unperturbed ground state |$o)- The fermion spinors are normalized 
according to w\ p w x , pl = <5a,a'<5 p ,p' and they obey the free Dirac equation 

[a ■ p + 7 m] w x , P = £\, P W\ iP (29) 

The interaction hamiltonian is 

Him = 9 J 0(x)^(x)^(x)rf 3 X 

= -7= E -7^«W+kWA,p&t, p+ k&A,p (°k + a -k 
V ft p>kjV VMs 



(30) 



We shall show that, in the ultrarelativistic limit, the quantum states correspond- 
ing to soft fermionic excitations can be obtained by diagonalizing the total hamil- 
tonian H = H + Hi nt , where H is the free hamiltonian for electrons and photons, 
in the following basis of states 

h, P \%) &„ lP+k a£|$o), (31) 

where ca, p < //, e^p+k < A* an d, as sa id earlier, |$ ) is the unperturbed ground 
state with energy E . All the states (|3ll) have momentum —p. The first state 
is a hole with momentum —p. The second represents a state composed of a hole 
with momentum — p — k accompanied by a photon with momentum k. These are 
eigenstates of H with eigenvalues E — e\ iP and E — e^ iP+ k + oJk respectively. A 
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pictorial representation of these excitations is given in Fig. 5. The single particle 
energy eA, p is, as mentioned above, positive or negative depending on A, namely 
£a,p = e P for A = 1, 2 and e x<p = — e p for A = 3, 4. 

We call X Xp and Y up ]l , respectively, the amplitudes associated with the two 
types of states (|3~ID in the decomposition of the eigenstates |^_ p ) of H on the 
unperturbed basis 

I*-p> = E X AA P |$o) + E^pA,P+k4l^o}. (32) 



Note that in restricting ourselves to the states (pT|) we are performing a truncation of 
the Hilbert space. In particular, the states with one hole and one photon are coupled 
to states involving two photons, two holes and one electron above the Fermi sea. This 
truncation is in line with the high density limit and is legitimate since the neglected 
states would contribute in higher order in g to the properties of the single hole state 
that we are interested in. The eigenvalue problem (H + H int )\^- P ) = — u\^- p ) 
then takes the form 

{uj - e x . P )X x , p -gY, J-n W,, p+k «;A,pK,p,k = 

~ y/2u Q ^ ^A,p^,p+k^A, P + {uj - e^p+k + Wfc)K )P) k = (33) 

where —uj measures the energy of a state with one particle less than |$o) with respect 
to E (see the discussion at the beginning of section 3). In the first equation, the 
sum runs over all single particle states such that e^p+k < f-i. 

In order to calculate the eigenvalues u, we use the second equation to eliminate 
the amplitudes ^, p ,k- One gets, after a straightforward calculation, 

(u - e A , p )X AiP - Xy, p w X ',pM(p, uj)w x , P = 0, (34) 

A' 

where M(p, u) is given by 

2 r d 3 k 1 f ^-p+k/p+k ^-p+k 



M(p,u) = g 2 P+KJP ! K + - P+ \ 70- (35) 

Here e p+ k = |ei/,p+k| an d Ap, introduced in Eq.(|2]) projects onto states with posi- 
tive or negative energy, i.e. A+ = J2u=i,2 w u, v w ] u p , and A p = J2u=3,4 w ^p w I,p ( see 
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Eq. (0)). The occupation factor fp + ^ is 1 for a particle in the Fermi sea and 
otherwise. The expression fl35|) coincides with the sum of the terms (ii) and (iii) in 
the expression ( |2"0"D of the mass operator S(p, uj). The term (i) in ( [2"0"D accounts, for 
small p and uj, for modifications of the Fermi discussed in the beginning; of 

section 3. Such modifications are negligible in the ultrarelativistic limit, and indeed 
they are not taken into account by the eigenvalue problem fl3~3"|). We shall come back 
to these at the end of this section. The last term of ([35]), which involves positron 
intermediate states, is also negligible in the ultrarelativistic limit, and we shall con- 
sequently ignore it in the following. By multiplying Eq. (|34f) by W\ tP and summing 
over A, one easily gets 

7 oG- 1 (p,^)£X* ipW A,p = 0. (36) 

A 

The values of uj for which Eqs. ( |3"3"| ) are satisfied are therefore the zeros of G" 1 , 
that is, as announced earlier, they coincide with the poles of the Green function 
calculated here in the ultrarelativistic limit. 

Using results established in section 2, we can write Eq. (|36|) as follows 

(w - a)x = [(1 + b)at ■ p + m^} x, (37) 

where x — J2\Xl p w\ tP . The eigenvalues are given by Eq. (|T2"D. Let us recall that 
there are two pairs of solutions which depend on the parameter e = ±1, with the 
normal solutions such that uje > 0. Besides, each solution possesses a two-fold spin 
degeneracy since the spin operator S = 75Q: commutes with both a ■ p and 70. 

The spinors W\ :P are eigenstates of the free Dirac hamiltonian (see Eq. (^)). In 
either case, p = or m = 0, one can choose the solution of Eq. (^) to be one of these, 
i.e. choose x — w x,p- Namely, when m = 0, we can choose x to be an eigenstate 
of a. ■ p with eigenvalue e. When p = 0, we can choose x to be an eigenstate of 70 
with eigenvalue e. In the general case where neither m nor p vanish, the presence 
of the factor (1 + 6) in Eq. ( |3"TD induces a mixing which depends on uj and which 
needs to be calculated explicitly for each eigenvalue. This will be done later, and 
we shall consider first these particular cases where m = or p = 0. Then, for each 
eigenvalue, only one amplitude, say X\ tP , is non vanishing. Setting X\ tP = 1, i.e. 
deferring the normalisation of the states to the end of the calculation, we obtain the 
following simple expression for the amplitudes l^ iPi k 

v - 9 ^,p^,k ( , 

v 2kQ t*J p cos 9 
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where we have used the fact that k > m to simplify the denominator, and 9 is 
the angle between k and p. 

In order to calculate explicitly the overlap W\ tP w u ^, we choose the chiral repre- 
sentation of the 7 matrices in which the spinors take the form 



M A,p = 



(e p + m + cr • p)(/> Ap 
2yje p (m + e p ) V (" e P - m + a ■ p)0 Ap 



2^e p {m + £p ) V (e P + m + • p)0 Ap J 

where p = p/p. These spinors are eigenstates of helicity with helicity A if the two 
components spinors are, i.e. if (cr • p)<p\ p = A0 Ap . In the limit m = they become 
eigenstates of 75, with eigenvalues A for u\ and —A for v\. In the limit p = they 
are eigenstates of 70 with eigenvalues +1 for u and —1 for v. Note that in both cases 
the spinor -u Ap is associated to e — 1 while the spinor v X p is associated to e — — 1. 
Since we now distinguish explicitly positive and negative energy solutions, the index 
A takes only two values, i.e. 1*1,-1 = ^1,2, v i-i — w 3,4- Accordingly, we change the 
decomposition of the fermion field into 

^(x) = 77^E(e~ ip ^A,p&A,p + e^t; A , P 4,-p)' ( 40 ) 



A,p 



where d\ _ p denotes the creation operator for a positron with helicity A and momen- 
tum —p. Also, from now on, |$ ) denotes the Fermi sea (i.e. excluding the Dirac 
sea). 

Let us first consider the case m — 0. Then we have 

UA lP u„,k = 4>\ p 4> v k^-\ UA )P w V) k = -A0l p 0^^, A . (41) 

Only spinors with opposite chirality contribute in the coupling to a scalar field. 
From these relations, one immediately deduces the form of the eigenstates 



i*«_ p)± , x = ^ 1 4,_ P + >e 1 '*->• (42 » 
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where |^" p ) (resp. |^!ip)) is the eigenstate corresponding to e = +1 (resp. e = 
— 1), the subscript + (resp. — ) refers to the normal (resp. abnormal) branch, 
and z± is a normalization constant. The norm of the states can be calculated 
directly from Eq. (f42|) by replacing the sum over k by an integral, as in Eq. fl35|) . 
and noting that <P\p4> X k = cos (6 l /2) and <j)\p(p_ x j t = sin(#/2). One then obtains 
z± = (u± — p 2 )/2M 2 , which coincides with the residues given by Eq . (|18D . It may 
also be verified by an explicit calculation, using the results of section 2, that all the 
states are orthogonal. 

One sees on Eqs. ( |i2"D that, as a result of the interactions, the "bare" hole &A,p|^o) 
and the "bare" positron d\ _ p |$o) couple independently to a a coherent superposition 
of hole-photon states. There is no direct mixing between the hole and the positron 
states. Such a mixing does in fact occur when m/0 and p ^ 0, as we shall see later, 
but even then, this is not a dominant feature. The basic mechanism at work is best 
understood by concentrating on what happens to one of the single particle states, 
the hole state for example, as the interaction is switched on. This is illustrated in 
Fig. 7 where, on the left hand side is drawn the dispersion relation u(p) for the hole 
state in absence of interaction, while in the right hand side is the split dispersion 
relation uj±(p) which results when the interaction is taken into account. The shaded 
area represents the phase space occupied by the continuum of hole-photon states to 
which the hole couples. We have represented only the dominant part of the density 
of states. As we have seen in section 3.1, the density of state is non vanishing in 
the region p < uo < e p , but it does vanish as ui — > e p . As it is clear from this figure, 
this coupling to the continuum shifts some of the single particle strength from the 
positive energy, where it normally belongs, to the negative energies where it appears 
as an "abnormal" branch. In the limit where m,p <C //, the strength becomes 
equally distributed among the two branches. The branches corresponding to the 
positrons are symmetrical to those just described for the hole. The interpretation 
of the two branches with positive u is then clear. The normal branch is that part 
of the electron strength which is pushed up by the interaction with the continuum; 
the abnormal branch is that part of the positron strength which is pushed up to 
positive energies by the same interactions. We note that the energy shift, being of 
order g[i rather than g 2 /i, is non perturbative. Most of the electrons in the Fermi sea 
contribute to this phenomenon, and indeed the structure of the states just described 
has many features of collective excitations in many body systems. 

We now discuss some properties of the spectrum in the light of these remarks. 
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When m = and p — > 0, — > M (see Eq. (p~5|) ) and z± — > 1/2 (see Eq. (0)) and 
all four states have a similar structure, i.e. the single particle strength is equally 
distributed among the four of them. When p increases, the abnormal branch quickly 
becomes uo_ ~ —p with a residue Z- — > 0, while the normal branch goes over to the 
usual fermion dispersion relation with z + — > 1. At the same time, the relative 
weight of the hole-photon component, (1 — z±)/z± decreases for the normal branch, 
but increases for the abnormal one. As revealed by a simple analysis, when p ^ M, 
only the states at the edges of the hole-photon continuum contribute to the abnormal 
branches. These states are of the form ^p+k^kl^o)? with p and k antiparallel when 
uj ~ —p, and with p and k parallel when u ~ p. 

The two states fl42|) , i.e. the dressed hole and the dressed positron, have the 
same fermionic charge —1 with respect to the unperturbed ground state. This may 
be verified explicitly by taking the expectation values of the operator Q = J : x/j'ip : 
<i 3 x = J2\, P (^a,p^a,p — ^a,p^a, p ) in either state. One may write this expectation value 



as 



Q± = ~(z± + (1 - z ± )) (43) 

where the first contribution is the "direct" contribution, i.e. that of the bare hole 
(or the bare positron), while the second contribution is the "induced" part, that 
is the charge carried by the hole-photon component. When p grows, z + — > 1 and 
z_ — > 0, as already mentioned. Thus, as p grows, on the normal branch the direct 
component grows, while on the abnormal branch it is the induced component which 
grows. Thus, when p ^ M, the single particle state completely decouples from the 
abnormal branch. 

One gets another illustration of this phenomenon by calculating the expectation 
value of the current operator I = / ip^oapd 3 *. = J2x, P P i^x^Kv ~ d\ p d\^. In a bare 
hole state &A,p|^o) this is I = — p, opposite to that of the bare positron with the 
same momentum, i.e. ^_ p - This difference may be understood by noting that, 
viewed as a hole in the Dirac sea, the positron is associated to a single particle 
state whose energy, — e p decreases with increasing momentum p, whereas the single 
particle state corresponding to the hole in the Fermi sea has an energy +e p which 
increases with increasing momentum. This difference translates into a difference in 
the signs of the velocities de p /dp of the two excitations. Taking the expectation 
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value of the current operator in the states(^) one obtains 




dcos# 



cos 9(1 — cos 9) 



{uj± — pcos9) 2 



(44) 



where we have neglected p compared to k in evaluating the second terms. It can 
be verified that these expressions coincide with the velocity of the excitations as 
deduced from the dispersion relation uj±(p). For example, for the hole excitation, 
we have 



with z± = 1/ (1 — <9£ u (co>, p)/duj), S u = a + bp denotes the mass operator in the 
eigenspace corresponding to e = 1, and I± = —e(duj±/dp). By carrying out the 
integrals in (f44"|), and using the results of section 2, one can rewrite the currents as 
follows 



As we did earlier in the case of the charge, we can interpret Eq . fl46|) as the sum of 
a "direct" contribution, that carried by the bare single particle, i.e. — epz±, and an 
"induced" contribution coming from the hole-photon component of the states (f42| ) . 
The hole photon part of the quantum states gives a contribution to the current 
which is opposite to that of the current carried by the bare particle. When p — > 0, 
I± = — ep/3. Thus, at very low momentum, the sign of the current is the same as 
for the corresponding bare particle: the negative group velocity on the abnormal 
branch with uj > simply reflects the fact that the corresponding bare particle is 
a positron state dressed by the medium, in spite of the fact that it has positive uj. 
When p becomes large, the current goes to 1 on both positive uj branches. However 
this is achieved in different ways on the normal and the abnormal branch. On the 
normal branch, the residue z + — > 1 and, at large p the current is carried entirely by 
the bare hole. On the abnormal branch, z_ — > 0, and the current is no longer carried 
by the positron component, but only by the hole-photon component. This explains 
the minimum (for uj > 0) in the abnormal branch of the dispersion relation. Indeed, 
as we have seen, at small p the current is mostly that of the bare positron. As 
p increases, the weight of the hole-photon component increases until it completely 
dominates. At this point, the current is that of the hole in the hole-photon state, 




(45) 



ep 2z± + — (l-2z±) 

L p 



(46) 
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and we have seen that the states with 00 ~ p correspond to holes with momentum 
parallel to p, and therefore to current opposite to that of the positron. 

We now turn to the case of a massive fermion at rest (p = 0). Then 



Ti 



V2 



At* 



(47) 



The eigenvectors are therefore 

to±,> = 




9 



E 



u±\/2 yk \/2uj k Q 
9 y -^Px<P„k 




(48) 



0. However, 
2 



The normalisation constant z± can be calculated directly, as for m 
both helicity states now contribute to the sum and one must use J2i. 
1. This gives immediately z^. 1 = 1 + M 2 /u±, in agreement with Eq.flTDp. The 
orthogonality of states with different u can also be checked explicitly. One could also 
repeat here the discussion given above in the case m = concerning the decoupling 
of the single particle excitation which occurs when m ^ M. 



A remarkable fact in Eq. fl48[) is that u factorizes out of the hole-photon com- 
ponent, so that, up to a trivial multiplication factor l/u±, the same coherent su- 
perposition of hole-photon states contributes to both the normal and the abnormal 
branch. Thus for example, we may write the eigenstate corresponding to the normal 
branch uj + in the form 



TO 



z + b Xfi \%) + Jl -*+!$!) 



(49) 



where |<3>i) is the collective state. The eigenstate corresponding to the abnormal 
branch uo_ is also a linear combination of 6a,o|^*o) an d \^xj- Since uj + 7^ the 
eigenstates are orthogonal to each other and we may write, up to a phase, 



TO- 



(50) 



and thus z_ = 1 — z + . The explicit construction of the eigenstates thus explains 
naturally why the residues add up to unity when p = (see section 2). On the 
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other hand, for a massless fermion with p ^ 0, Eq.(|42|) shows that the hole-photon 
component is different for the two branches and the previous argument cannot be 
used. Indeed, the sum of the residues is smaller than unity in this case, some of the 
specral weight being displaced in the region — p < uj < p. 

Let us finally consider the general case when m and p are both different from 
zero. Then the solutions x of Eq.fl3~7|) are linear superpositions of the solutions w\ jP 
of the free Dirac hamiltonian. For instance, a solution with e = +1 and helicity A 
can be written as a superposition of u\ tP and V\ tP defined in Eq.(p^): 



X 



1 - r± m a ,p + e v y/r± v x , P 



(51) 



where r + (resp. r_) measures the mixing between the hole state ^a.pI^o) an d the 
positron state <i A _ p |$ ) on the normal (resp. abnormal) branch. For a state with 
e = —1, u and v must be interchanged in the latter equation. Thus the eigenstates 
are of the form 



l*- P > 



Z± V 1 - r ± & A,p|$ ) + 4,-pl $ o) + V 1 - Z±\®l) 



(52) 



if e = +1, and b and <$ must be interchanged if e — — 1. In Eq.(^), z± is the residue 
at the quasiparticle pole and is a linear superposition of states fe^p+k^kl^o) as 
in Eqs.([42"D and (^). From Eq . (|52|) we deduce that the probability that a hole state 
^A.pl^o) be found in the eigenstate |\l/" p ) ±jA (resp. |^ r !ip}±,A) is ^±(1 — r ±) (resp. 
z±r ± ). 

The mixing coefficients r± can be calculated easily by noting that the matrices 
7° and a. ■ p appearing in Eqs. ( f2~9|) and (|37|) obey the same commutation relations 
as two arbitrary Pauli matrices, say o\ and o"2. On may then associate a solution to 
Eq. ( p9|) to a spinor pointing in a direction making an angle 7 with the 1-axis, with 
tan7 = p/m, and the corresponding solution to Eq. (37) to a spinor pointing in a 
direction making an angle 7' with the 1-axis, with tan7' = (1 + b)p/m. The overlap 



between the two spinors is then 



p 



\J\ — r± = cos((7 — 7O/2) or, explicitly, 



1 + (1 + b)p 2 /m 2 



\Jl + p 2 /m 2 ^Jl + (1 + b) 2 p 2 /m 2 



(53) 



It depends on uj through b. For the normal branch (resp. the abnormal branch), 
b varies from to —1/3 (resp. from —00 to —1/3) when the chemical potential 



26 



increases from to +00. Thus the bounds on the mixing coefficient r± are 



< r 1 < 



1 + (2/3)p 2 /m 2 



^1 +p 2 /m 2 - v /l + (4/9)p 2 /m 2 



< r_ < 



p/m 



Jl +p 2 /m 2 



.(54) 



The maximum value of r + for the normal branch is reached when p/m = y 3/2 and 

infinite density: r™ ax = 1/2 — y/E/5 ~ 10~ 2 , which is rather low. On the other 
hand, for the abnormal branch the mixing can exceed 1/2 if the density is not too 
high; however, the residue z_ of the abnormal branch is then very small, so that 
the product z_r_, which is the probability that a hole state 6a, p |$o) be found in the 
state \ty v )_ t \, remains small, at most of the order of 10~ 2 . 

This mixing between holes and positrons can be easily understood in the ultra- 
high density limit M 3> m,p. In this limit Eq.(|T2"D gives uj — > ±M for both values 
of £ (see Fig. la). Thus there is a degeneracy at uj = M (resp. — M) between the 
branch uj + (normal hole branch) and the branch — u;_ (abnormal positron branch). 
This degeneracy is lifted by terms proportional to p and m, which act like a pertur- 
bation in the eigenvalue equation (|3~7"D and therefore induce mixing between the two 
branches. Writing Eq.(|37|) in the form (H (uj) + V)x = ujx with H (uj) = a ~ M 2 /u 
and V = (2/3)a ■ p + 7°m (we have used the fact that b — > —1/3 for ^> p), and 
expanding around u = M one gets 



'dH 



M) + V 



LU = M 



X = (u- M) X - 



(55) 



Since (dHo/du) w =M - 
one gets immediately the eigenvalues uj 



1, the problem reduces to diagonalizing V . In particular, 

M ± y/p*/9 + m 2 /A. Note that r+ = r_ 
in this limit, and that the mixing vanishes if either m = or p = 0, as expected. 



The calculation of the eigenstates that we have performed in this section, which 
is equivalent to solving Dyson's equation for the poles of the Green function in the 
ultrarelativistic limit, assumes that the modifications of the ground state can be 
ignored. More generally, the solution of Dyson's equation for the full one loop self 
energy is equivalent to the following eigenvalue problem 



eA,p)^A, P 



9 



9 



kv 



-kW A , p X A . p + (u 



w 



-k^A,p^,p,k — 
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~ — H ^, P +kWA,p^A,p + (o> - e^,p+k - ^fc)^, P ,k = (56) 



where, in the first equation the sum runs only over the states such that e p+ k < \x for 
the Y amplitudes, and such that e p+ k > A* f° r the Z amplitudes. That this system of 
equations is equivalent to Dyson's equation can be easily verified by eliminating the 
amplitudes Yi,,p,k and Z v v ^ with the help of the last two equations, and comparing 
the resulting equation for X\ v with Dyson's equation. This comparison allows us 
to identify the mass operator, and one recovers indeed Eq. (pop . The equations (|56|) 
take into account the modifications of the Fermi sea due to interactions, at this order 
of perturbation theory. In particular, because of the interactions, single particle 
levels are not fully occupied below the Fermi level, and there are non vanishing 
amplitudes of the form 

2*,p,k = (^. P ,k|6i,p+ k ak|*o), (57) 

where |^ ) an d (^, P ,k| denote true eigenstates, that is including perturbative cor- 
rections. As for the amplitudes X and Y they are given as before by 

X u ,k = (^u,pK P \^ ) Y v>Pjk = (* V)P ,k|&»,,p+k4l*o}, (58) 

but now, the state \^>q) contains perturbative corrections, i.e. it differs from the 
unperturbed Fermi sea |$o)- It is easily seen on Eq. (|56|) that in the ultrarelativistic 
limit, and for soft modes, the amplitudes become negligible as compared to 

the y^ p ,k ones. 



5 Discussion 



We have, in the main text, described a model plasma composed of Dirac fermions 
interacting with a scalar field. Most of our results extend however to a dense system 
of electrons (see Appendix C), or to dense quark matter. In such high density 
systems, particle masses are negligible in a first approximation, and the only energy 
scale is provided by the chemical potential, [i. Most particles in such a system have 
momenta of order /i. Long wavelength excitations develop on the scale g/j,, and have 
a collective character. Bosonic modes, not discussed in this paper, correspond to 
the usual plasmons involving particle-hole excitations localized in the vicinity of the 
Fermi surface. We have concentrated here on fermionic modes, which can be excited 
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by processes adding small momentum positrons or removing electrons at the bottom 
of the Fermi sea. 



The physical picture which emerges from our calculation is quite in line with 
the work of Ref . |13| where the equations of motions describing the long wavelength 



excitations of an ultrarelativistic plasma at high temperature are reduced to a set of 
coupled mean field and kinetic equations incorporating the dominant medium effects 
as "polarization" phenomena. In either case of large temperature or large chemical 
potential, the polarization is entirely due to the hard particles with momenta of 
order T or fi. In the zero temperature case, the dominant process is that in which 
a hard electron annihilates with the soft hole, or the soft positron, and turns into 
a hard photon. Thus, the quantum numbers of the hard particles fluctuate, with a 
typical period of order This is accompanied by a corresponding oscillation in 

the number of soft particles, which is interpreted as an oscillation of the average soft 
fermionic field in Ref. The present study allows us to give a simple illustration 
of this in terms of the quantum states that we have constructed. Consider indeed a 
perturbation which consists in adding soft positrons to, or creating soft holes in, the 
unperturbed Fermi sea at time t = 0, say. The system will then evolve into a state 
|$(£)) which will contain components on the states \^ £ )± constructed in section 4: 

\$(t)) = \%) + u(tW)± (59) 

where a(t) is a small time dependent amplitude. The expectation value of the 
fermionic field in such a state is proportional to the amplitudes X\ p introduced in 
section 4, while the amplitudes Yx^ ^ contribute to the "induced source" for the 
fermion field. 

The explicit construction of quasiparticle states in section 4 clearly shows how 
new modes emerge: the free electron hole is coupled to a large number of states 
which do not interact together and have an energy comparable to that of the un- 
perturbed hole excitation. To outline this structure, let us consider a simplified 
situation where a quantum state |0) with energy E (corresponding to e p in the 
previous section) is coupled with equal strength V to a large number N of states 
\i), uniformly distributed throughout the energy interval [— AE/2, AE/2]. In the 
calculation performed in section 4, AE = 2p, N is of order Qfi 3 in a volume Q and 
V is of order g/^/Vtji, i.e. Vy/~N ~ g\i. The equation determining the modes uo is 



NV 2 

u-E = -—log 



uo + AE/2 



uo- AE/2 



(60) 
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A substantial energy shift and/or a new quasiparticle will occur only if V^/N is at 
least of the same order of magnitude as both Eq (which is the difference between 
the energy of |0) and that of the centroid of the states |z) to which |0) couples) and 
AE (the energy interval over which the states \i) are spread). For the calculation of 
section 4, this is equivalent to saying that g\i must be at least as large as m and p. 

One can then wonder whether a similar situation could occur in other systems. 
In fact, a structure similar to that studied here has been found in the study of the 
electron gas at intermediate densitites[[14|, [TH] . There, the coupling of an electron 
hole to collective plasmons gives rise to a peak in the imaginary part of the hole 
self-energy which is sufficiently strong to produce a new quasiparticle, dubbed the 
"plasmaron" . However, in the case of the non relativistic electron gas, the validity 
of the perturbative approach is less clear, and it was indeed found that vertex 
corrections could very well suppress the structure J0J. 
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APPENDIX A : The vacuum part of the self— energy 

At zero temperature and zero chemical potential, the mass operator So(p, a?) 
may be written as 



S (p,w) = a (p 2 ,p )j + 6 (p 2 )7 ■ P + co{p 2 



(A.1) 



where, because of Lorentz invariance, the functions ao/po, bo and Co depend only on 
the 4-momentum squared, p 2 = pi — p 2 , and b (p 2 ) = — a (p 2 ,Po)/Po- We shall limit 
the discussion to the case p = (the result for p ^0 may be obtained from Lorentz 
invariance). The expressions for and Cq may be read off Eqs. @, and they involve 
ultraviolet diverging integrals over momenta. Adding the usual counterterms of mass 
and wave function renormalization, we write the renormalized inverse propagator as 



= G + T,q(u) — 5m 
= (-lu + a (uj)) 7° + (m + Cq(u) - 5m) 
The counterterms are fixed by the on-shell renormalization conditions 



(A.2) 



5m 
Z^ 1 - 1 



a (m) + c (ra) 



1 - 



da 



did 



dc 
dcu 



(A.3) 



In doing the calculations of Z 2 , it is necessary to give the photon a small mass A 
to avoid the infrared divergence. The Green function can finally be written in the 
form 



-lo + a R {uj) + m + c R (u), 



(A.4) 



where a R (u) and Cr(u) are the renormalized contributions of the vacuum to the 
mass operator, which must be added to a\(u) and c\(u) in the finite density case: 



a R Uo) 



cr(uj 



16vr 2 
g 2 m 

167T 2 



<jj- 



m 



2uo 



1 + 



m 



In 



ijj- 



1 - 



- u 5 + 4 In 



A 



/ m 2 \ , 










+ (5 + 41n- j 






y my 



(A.5) 



and we have kept only the dominant terms in X/m as A — » 0. 
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Note that the renormalization does not affect the imaginary part which is finite. 
Note also that in calculating the frequencies of the modes at high density, one can 
ignore a R and cr. Indeed, when uj ~ gfi, their contribution is of order g^fi while that 
of Si is of order gfi . Thus, in particular, the term in A which is infrared divergent 
and, in gauge theories, gauge dependent, does not contribute in the high density 
limit. Presumably it would disappear in a consistent calculation of quasiparticle 
energies at order g 3 fi, i.e. beyond one-loop . 



32 



APPENDIX B : Exact expressions for the one-loop mass operator at 

T = 

In this appendix we give the exact expression of the mass operator (EDI), valid 
when either p or m vanishes. 

If p = 0, the finite density contribution to the mass operator can be written 



S 1 (p = 0,w) = a 1 (cu)7° + c 1 (cj) 
The real parts of a\ and c\ are 

'VP 



(B.l) 



Re [ai(o;)] 



Re [ci(cj)] 



g 2 m 
16tt 2 



g 2 m 
8^ 



m 



m 



l ( P 



1 + — 7 + - + — o cosh — 



+ 



2m 



1 - 



m 



r 



In 



u — fJ, 



m 
k 



uj — n + \J [i 1 — m 2 



hp m 2 , / u 

— + — cosh 1 ( 



1 - 



m 



In 



UJ — fi+ y/ fi 2 — m 2 



(B.2) 



with kp = kp the Fermi momentum. The imaginary parts are (including the vacuum 
contribution) 



Im [di(co')] 
Im [cx{u)\ 



-7T 



2 4 4 

<7 m — 



-7T 



167T 2 2cj 3 

2 2 2 

g m m — uj 

167T 2 UJ 2 



(B.3) 



if 



< — m or fx — kp < uj < m 



(B.4) 



and otherwise. Note that this equation agrees with Eq . (P2|) when we set p = in 
the latter. We note the relation Im(a + c) /Im(a — c) = (uj + m) 2 / (uj — m) 2 , identical 
to that found at high T and // = |. The imaginary part is independent of /i, 
except for the domain where it is non vanishing. In the high density limit, this 
domain is the interval m 2 / (2/x) < uj < m; the limit m — > is thus singular. 
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If m = 0, hp = Hi the finite density contribution to the mass operator can be 
written 



£i(p, u) = ai(p, uj)<y° + 61 (p, uj)j ■ p 



(B.5) 



The functions a\ and b± can be calculated analytically by choosing (p — k) as inte- 
gration variable in Eq. (0). The real parts are: 



Re [ai(p,w)] 



32vr 2 p 



p^+^-^-^^ + ^+^ln 



+ p 



(2/i - a; + p) (2/i + - p) 



In 



2fi — lu — p 
to — p 



2/i — uj + p 



Re [fei(p, a;)] 



167T 2 p 2 



2 , 2/i-uj-p 

2/i — u; + p 
4p 



a; 



cu 2 -p 2 - -(2/i + cu + p) 



a; 



uj — p (2/i + a; — p) 



In 
In 



w + p 



2/i — uj — p 
uj — p 



2/i — uj + p J 
(B.6) 



As to the imaginary parts, several cases must be distinguished: 

,2 



in 



II : 

IV : 
la : 



lb : 



Im [a(p, uj] 
Im [b(p, uj] 
Im [a(p, uj] 

Im [5(p, uj] 
Im [a(p, a;' 
Im [a(p, uj] 

Im [5(p, a;' 

Im [a(p, uj] 
Im [5(p, uj] 



9 

9 2 
- n 32^ 

g 2 (2/j,-uj- p) (2/i + u + p) 
32tt 2 4p" 
g 2 2/i — uj — p f , a; 



-7T 



-7T- 



-p - -(2/i + a; +p) 



32vr 2 2p 3 
Im [5(p, a;)] = 

g 2 (2fj,-uj-p)(2fj l + u + p) 
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32n 2 4p 
g 2 2/2-uj-p 



32tt 2 2p 3 
9 2 

9 2 
"32^ 



UJ 



uj p -(2/i + a; +p) 



(B.7) 
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where the various regions are indicated on Fig. 10 (III : uj < —p;II ■ —p < 00 < 
min(p, 2/i — p); IV : min(p, 2/i — p) < uj < max(p, 2\i — p); la : max(p, 2fi — p) < uj < 
2/i+ p; lb : u> > 2/i+p). These expressions include the contribution of the vacuum to 
the imaginary parts. One checks easily that these equations give the same result in 
the limit p — >• as Eqs.( [B.2| ) and ( |B.3| ) when m — > 0. Furthermore, if \x 3> ou,p, one 
recovers Eqs.(@) and Eq. (p70|) . Note that the imaginary part is everywhere negative. 
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APPENDIX C : QED 

The results obtained in sections 2 and 4 in the case where fermions are coupled 
to a massless scalar field can be easily generalized to an electromagnetic coupling. 
The corresponding modifications are summarized in this appendix. 

In the Coulomb gauge, the electromagnetic field is decomposed in 

1 ' * 1 ( —ik-x M , ik-x _t 



#(.*) = ^ £ p E ^ + ^' «S«W . (c.i 

where p = 1, 2 is a polarization index and k are normalized polarization vectors 
e pk = 0' e p,k • k = and e P)k ■ e p ' ik = <5 PiP /. The interaction hamiltonian is now 



ffmi = g J 1p(x)fi(x)lp(x)d X 
V " p,k,A,i/,p V fc 

(C.2) 



and the result of the one loop calculation may be obtained, in the ultrarelativistic 
limit, by diagonalizing H in the basis 



9 A,p 



$o> ^, p+k aL|$o). (C.3) 



The corresponding amplitudes will be denoted by X A p and ^A,p,p,k- The eigenvalue 
equation takes the form ([34]) with 

* A+ ^,0J . f .. </* A- o,0 i . 

(C.4) 



-x 2 f d 3 k 1 f ^, k A+ +k 7Vp,k/ P +k ^, k A p+k 7Vp,k 

M{p,u)=g / + _' ■ P 

J (27T) d 2a> fc I -a; - e p+k + u k -u + e p+k + u k 

In the ultrarelativistic limit (k 3> m,p), the term proportional to m in A ± (p + k) 
can be neglected, and p + k replaced by k. Then, since e Pik • k = 0, / Pjk = —7 • e P)k 
anticommutes with 7 ■ k and thus with A k 7°. Then, using Dp^p, k /p,k = —2, one 
gets 

M(p, d>)coulomb = 2M(p, u;) scalar . (C.5) 

Thus, the only difference with the Yukawa coupling is a global factor of two in the 
mass operator, which corresponds to the fact that there are two spin states for the 
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bosons instead of one. The same result holds more generally for the mass operator 
E in the ultrarelativistic limit, even at finite T, so that the results of section 2 are 
also valid in the case of QED, if the definition of M in Eq.(0) is replaced by 



M = ^yV 2 + k 2 T 2 . (C.6) 

Let us now construct the eigenstates explicitly. The difference with the scalar 
field lies in that the spinor overlap w u ,p+kWx )P in Eq. (|33|) must be replaced by 
M^p+k^k^A.p- If rn, = 0, 

^A,p^, k ^,k = -Mxptp^ ■ &<i) v k$v,\ v X;P £* tk u Uik = -0lpe Pik ■ (T(j) u% 8 v - X . (C.7) 

The electromagnetic interaction, unlike the scalar coupling, conserves chirality. The 
eigenstates are therefore 

*- p /±,A - V Z ± °A, P + 9 L^~rK=R a — °A,p+ka„k \®o) 

p \ V 2u} k Q w± - p cos6» p ' J 

U r,v \ \^ 1 ~^Ap e p,k • °^-Xk h f 

p v i , p ^ v/2u; fc fi w± - p cosfl p ' 



In order to calculate the normalization, the electric current, etc, one must evaluate 
the sum over polarization states J2 P 0Ap e p,k ' <T 4 > ±xk = 1 ~F cos ®- These expressions 
are similar to those obtained in the case of a Yukawa coupling, except for a factor 
of 2 which is absorbed in the definition of M. So the results are finally the same. 

In the case p = 0, a similar calculation gives 

„,<< - — | , 9 sr ~ v ' °hk h \\ iris \ 

^±v ^ v^fc" y 

I*-p) ±j a = ^ | <o + 77^ E ~^£i k n* vi K *<*) |(&o) - (a9) 
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FIGURE CAPTIONS 



Fig.l: 

The spectrum of fermionic excitation in an ultrarelativistic plasma for two cases, 
M = 10m (Fig. la) and M = m (Fig. lb). The full lines correspond to e = 1, the 
dashed lines to e = —1. The lower part of the figures displays the residues at the 
quasiparticle poles: z + is the residue of the normal branches and Z- is the residue 
of the abnormal branches. 

Fig.2: 

Graphs corresponding to the three contributions labelled (i),(ii) and (iii) in the 
expression (|20|) of the mass operator. 

Fig.3: 

Energy domains where the imaginary part of S, given by Eq. (p0|) , is non vanishing. 
The three regions (I), (II) and (III) correspond respectively to the electron-photon 
continuum, the hole-photon continuum and the positron-photon continuum, and are 
the regions in which, respectively, the contributions (i),(ii) and (iii) to ImS are non 
vanishing. 

Fig.4: 

Diagrams contributing to the second order correction to the energy of a hole (i) with 
momentum p < kp , and a particle (ii) with momentum p > kp, and reflecting mod- 
ifications of the Fermi sea. The diagram (0) represents the second order correction 
to the energy of the Fermi sea, i.e. AEq (see text). It is the sum over the fermion 
momenta in this diagram which is restricted by the presence of the extra particle 
(in (ii)) or the extra hole (in (i)). 

Fig.5: 

Schematic representation of some important excitations of the ultrarelativistic plas- 
ma at zero temperature, (a) A hole at the bottom of the Fermi sea. (b) A positron 
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with p = 0, i.e. a hole at the top of the Dirac sea. (c) A hole left by a "hard" 
electron with momentum p ^ kp. This hole is accompanied by a hard photon (not 
drawn) . 

Fig.6: 

a) solid lines: real part and imaginary part of the mass operator for an eigenstate of 
7° with the eigenvalue e = 1. The calculations are done using Eq. ( |B.2| ), Eq. ( p.3|) 
and Eq. ([Bl]) (the vacuum contribution to the real part is thus omitted). The 
imaginary part vanishes stricly on this graph for — m < uo < and for uo > m. 
It is small but not zero for uo < —m. The labels on the curves (5,10,20) are the 
corresponding values of /i/m. In this figure as well as in Figs. 4 and 5, g 2 / An = 1/137. 
Vertical dash-dotted lines: indicate the value u> = /i — kp where the mass operator is 
divergent. The other divergence, at uo = n + kp is outside the range of the abscissa. 
Dashed line: line uo — m, whose intersection with the solid line gives the elementary 
modes according to Eq. (|23|) . 

b) Same as a) for an eigenstate of 7 with the value e = — 1. The dashed line is now 
uo + m according to Eq. ( |23| ) . 

Fig.7: 

a) solid lines: real part and imaginary part of the mass operator for an eigenstate 
of a • p with the eigenvalue e = 1 . The calculations are done using Eq. ( |B.6|) and 
Eq. ( p.7|) . The imaginary part vanishes strictly on this graph for uo > p. It is 
small but non vanishing for uo < —p. The labels (10,20,40) on the curves are the 
corresponding values of fi/p. Dashed line: line uo — p, whose intersection with the 
solid line gives the elementary modes according to Eq. (|26|). 

b) Same as a) for an eigenstate of a • p with the value e — — 1. 
Fig.8: 

Same as Fig.7 with p very close to the Fermi surface {kp/p = 0.8, 1, 1.2). 
Fig.9: 
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The dispersion relation for the hole excitation. Left-hand side: no interaction. 
Right-hand side: with interaction. The shaded area represents the most important 
part of the continuum of hole-photon states to which the hole couples. 

Fig.10: 

The energy domains where the imaginary part of the mass operator is non vanishing, 
in the limit m = (see Eq.(B.7)). 
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